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Anomalously high generation efficiency of coherent higher 
field-harmonics in collisions between oppositely charged par- 
ticles in the field of femtosecond lasers is predicted. This is 
based on rigorous numerical solutions of a quantum kinetic 
equation for dense laser plasmas which overcomes limitations 
of previous investigations. 



Modern high intensity short-pulse lasers Q have made 
it possible to create dense correlated plasmas under lab- 
oratory conditions promising a large variety of applica- 
tions, e.g. Remarkable progress has been made in 
the field of time-resolved plasma diagnostics 0], inter- 
esting high-field phenomena have been observed, includ- 
ing production of multi-keV electrons and the generation 
of higher field harmonics (bremsstrahlung, BS) in neu- 
tral gases. On the other hand, despite early predictions 
||, harmonics generation in an ionized medium due to 
collisions of oppositely charged particles in the field has 
not been observed experimentally j6|. The reason is the 
low electron-ion collision frequency v el in a low density 
equilibrium plasma |Q. fn this Letter we demonstrate 
a drastic, more than 6 orders of magnitude, efficiency 
increase for dense nonequilibrium plasmas excited by a 
femtosecond laser pulse which should make this funda- 
mental effect observable experimentally. Moreover, we 
predict coherent harmonics pulses which may be signifi- 
cantly shorter than the pump pulse. 

A full theoretical description of dense laser plasmas 
requires the selfconsistent treatment of i) hydrodynamic 
mass and energy transport of electrons and heavy parti- 
cles, ii) collective (mean-field) phenomena and iii) elastic 
and inelastic carrier collisions in the presence of the field. 
In recent years there have been substantial advances in 
the theoretical description of the first two effects, see e.g. 
HHH and [^p|-p2|, respectively. However, these treat- 
ments, usually being based on Vlasov, Fokker-Planck or 
particle in cell simulations, either neglect collisions com- 
pletely or study them in a strongly simplified manner, 
fn particular, they neglect the field effect on the collision 
process. This may be justified for very strong fields where 
the electron oscillation velocity vq = eEo / (m e Q) is much 
larger than the thermal velocity Vth , and the collision fre- 



quencies v ee and v e i are much below the oscillation fre- 
quency f2 of the field. However, for field intensities below 
10 15 W/cm 2 ) and/or high plasma density (n > 10 20 cm~ 3 ) 
and low temperature, collisions are the dominant fac- 
tor in field-matter interaction at short times. At a later 
stage, due to hydrodynamic motion of electrons and ions, 
plasma heating and expansion, field effects will govern 
the plasma behavior, although collisional processes such 
as inverse bremsstrahlung, e.g. |I^J^,^] or impact ioniza- 
tion/recombination remain important. 

When this transition occurs and what physical pro- 
cesses take place at the early stage is not known un- 
til now - to answer these questions is the main goal of 
this paper. Our analysis reveals that during this ini- 
tial stage the plasma has a strongly anisotropic and non- 
Maxwellian electron distribution function which develops 
bunches of fast electrons. As a result of electron-electron 
collisions, eventually a monotonically decaying distribu- 
tion is approached, but this thermalization takes about 
trei ~ 50 . . . f 00 fs. Finally, the most interesting effect is 
that, under the strong nonequilibrium conditions at these 
short times, the generation of higher field harmonics in 
e-i collisions is strongly enhanced. 

The analysis of the initial stage of laser-matter inter- 
action requires a quantum kinetic equation which fully 
includes nonidcal plasma effects the importance of which 
has been noted frequently, e.g. [jsUI^l. Such an equation 
has recently been derived |Q and has the form 



d_ 

dt 



e a E(i)V k i/„(k a ,i) = ^7 ab (k a ,t), 



(1) 



where for the initial period spatial homogeneity may be 
assumed. The collision integrals are given by 

I ab (K,t) = 2j yV ^q)! 2 ^-^) 



X / dt COS 



(t-t)-^R ab (t,i) 



X {fa ft [I - fa] [I - fb] - fafb [1 - fa] [1 ~ fb}} | p (2) 

where a, b label electrons and ions, and we denoted e a b = 
( a + £fc, e a = pl/2m a , k ab = k a + k h , and q = k a - k a . 
V a b is the Fourier transform of the screened Coulomb 
potential V a b(q) — 47re a e(,/[g 2 + «(t) 2 ], [n(t) is the in- 
verse nonequilibrium screening length]. This allows to 
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avoid any cutoff (Coulomb logarithm). We underline that 
the collision integral I a b is exact (within the static weak- 
coupling limit). In particular, it permits a rigorous in- 
vestigation of the short-time physics and, for long times 
t 2n/ui p i [u)pi = (47rne 2 /r7i) 1 / 2 is the plasma frequency] 
and zero field (R a b = 0), it goes over to the familiar Lan- 
dau collision integral J dtcos... — > hir5(e a b — £ab)- More- 
over, this integral includes Pauli blocking (spin statistics) 
and it conserves the total energy (kinetic plus potential 
energy) which is important for dense nonideal plasmas. 

Most importantly, the integral I a t> includes the influ- 
ence of a strong field on the collision process of two par- 
ticles of species a and b exactly, as well as the effect of 
a finite collision duration r co n = t 2 — ti which is cru- 
cial for the results shown below. Indeed, during the time 
T co n each of the colliding particles gains a certain momen- 
tum Q a {ti,t2) in the field, and a pair of colliding (oppo- 
sitely charged) particles is displaced by the field a dis- 
tance R a b(ti,fo)- For a harmonic field, E(t) = E cosf2£, 



Qa{tl, t2 



(sin Qti — sin f^), 



(•3) 



Ra&(*i!*2) — ( — - r ) ?>§ 

x {cos flti — cos Qt 2 — fl (t 2 — 1 1) sin fit i} , (4) 

and these effects are growing with increasing field 
strength. The momentum shift enters the arguments 
of the distribution functions in the collision integral 
(0), (time-dependent intra-collisional field effect): f a = 
f [k a + Q a (t,t),t\; f a = f [k a + Q Q (t, t), t\, whereas R afc 
modifies the energy balance of the two-particle collisions 
in a monochromatic field. Thus the collision integral (|^) 
selfconsistently includes nonlinear field effects, such as 
harmonics generation, emission (absorption) of laser pho- 
tons during the two-particle scattering [(inverse) BS] and 
generalizes previous treatments of these phenomena, e.g. 
HP), to the case of arbitrary short times, nonequilibrium 
distributions and modifications by dense plasma effects. 

We expect the most interesting physics to occur in sit- 
uations where the effects of collisions and of the field are 
of the same order: (i) when the field frequency fl is in the 
range of the collision frequency v e i and the plasma fre- 
quency ujpi and, (ii) when the field induced particle veloc- 
ities vq are of the same order as the thermal velocity Vth- 
Thus, considering as an example a dense fully ionized hy- 
drogen plasma in a monochromatic optical field of inter- 
mediate wavelength A = 500nm (fl = 3.77fs _1 ), the fol- 
lowing parameter range is of interest for the simulations: 
densities 10 20 cm~ 3 < n < 10 24 cm~ 3 , corresponding to a 
plasma frequency range 0.56 fs _ < lu p i < 56 fs _ , and 
an initial plasma temperature of 20, 000-ftT which may in- 
crease up to almost 10 6 K within the first 50fs of the re- 
laxation. Correspondingly, for this case field strengths 
of interest are 10 7 V/cm > Eq > 10 9 V/cm. For stronger 
fields (/ > 10 15 W/cm~ 2 ) collisions will already be of mi- 
nor importance. 



We solve equations (|lj-|4|) by direct numerical integra- 
tion starting from a pre-excited electron-ion plasma in 
equilibrium, and the relaxation is computed over about 
50 field cycles. We underline that we solve equations (|l|- 
^) without any simplifying assumptions. In particular, 
momentum anisotropy of the electron distribution func- 
tion as well as time and momentum retardation are fully 
included. The ion distribution was found to remain a 
Maxwellian during the considered time interval which al- 
lowed us to simplify the integral I e i . To properly account 
for the time-dependence of the screened Coulomb poten- 
tial V(q), the inverse screening length n(t) is computed 
selfconsistently from the current nonequilibrium distri- 
bution functions. Solving the kinetic equation yields 
the time and momentum-dependent distribution function 
/ e (k, t) which allows us to compute all transport prop- 
erties, including the mean electron kinetic energy den- 
sity (E^ in )(t) = ^{v"l) and the electrical current density 
j(t) = X) a e a (v a ), [/ is normalized to the density, and 
(. . .) denotes averaging over /(*)]. 
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FIG. 1. Average electron kinetic energy increase versus 
density, with and without e-e scattering included. Field am- 
plitude and wavelength are E = 3 x 10 8 V/cm, A = 500nm, 
the initial plasma temperature is To = 20, 000K. 

Let us now come to the results. First recall that, in 
the collisionlcss case (I a b = 0), the electron distribution 
would oscillate as a whole according to sin Qt without 
changing shape. This yields no net current and no en- 
ergy increase (the energy oscillates with 2f2 around the 
ponderomotive energy U pon d = tou 2 /4). In contrast, in 
the presence of collisions, we observe a strong increase 
of the mean electron kinetic energy. After an initial 
transient of about 5 fs, (E^ in ) grows almost linearly, 
thereby oscillating with 2Q. This energy increase is due 
to field-induced electron acceleration during electron-ion 
collisions, cf. the term qR ei in the integral I e i : Eq. (|^). 
No net acceleration occurs in electron-electron collisions 
(R ee = 0). The energy growth rate depends on the ra- 
tio of field strength and density. This is shown in Fig. |l| 
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where the average slope of the energy curves is plotted 
for different densities. Interestingly, the most efficient 
plasma heating occurs around 10 22 cm -3 which is near 
the critical density n cr — 4.5 • 10 21 cm -3 where uj v i = Q. 
The maximum value of the power is about 1.5eV/fs which 
roughly corresponds to the absorption of one photon per 
electron and laser cycle. Notice also that, although e- 
e scattering does not contribute to the heating directly, 
this process may not be neglected, as this would lead to 
an overestimate of the power by a factor of three and a 
shift of the maximum to higher densities by one order of 
magnitude, cf. the triangles in Fig. [l]. 
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FIG. 2. Electrical current spectrum in a time-dependent 
electric field. The higher harmonics are the sole result of 
electron-ion collisions, (a) - time evolution of the spectrum 
for n = 10 22 cm -3 , from averaging over 5 laser periods (see 
text in fig.), (b) - spectrum (from averaging over the full 
calculation) for various densities. Other parameters are the 
same as in Fig. |lj 

Next, we consider the electrical current density. In con- 
trast to the collisionless case, where jo(i) = env sinfti, 
collisions between electrons and ions may give rise to 
higher field harmonics. Indeed, electrons oscillating with 
the field are disturbed by the Coulomb field of ions which 
yields a change of the electron velocity component in field 
direction (v > v th ) 



Av z (t) 



ay/2 



sin fit 



n(i + b 2 ) vi + 2b 2 



with a = e 1 , b — 




ro = 



cos 2ilt 
m e rrij 



(5) 



ft 2 ' 



where /i -1 = m ( r 1 +m,r 1 , and r± is the minimal electron- 
ion distance. One readily checks that in the spectrum 
of (||) appear all harmonics of the field. Formula (||) 
indicates that the spectrum will be particularly broad 
for close encounters (small r±, strong scattering). Fur- 
ther, the efficiency of harmonics generation decreases 
with growing field strength. On the other hand, for the 
case of a weak field [not covered by Eq. (g)], the har- 
monics yield increases with E (with the number of pho- 
tons), resulting in a maximum of the efficiency around 
^o — v th Quantitative estimates for the conversion 

efficiency rj n (intensity of the iV-th harmonics normal- 
ized to the intensity of the fundamental oscillation) can 
be given for the case of Maxwellian electrons. Then only 



odd harmonics exist and ru^ ~ (v e i/fl) 2 /N aN , where 
Ojv = a N(vo/vth) > (|. Direct evaluation of the 
collision integral I e i, Eq. (^), for a dense equilibrium 
plasma revealed that the highest value in equilibrium is 

rg Q ~ icr 7 (TJ. 

However, as we show now, the efficiency may be in- 
creased drastically under nonequilibrium conditions of 
femtosecond laser pulse excitation. Fig. || shows the 
results for the current spectrum obtained from the nu- 
merical solution of the quantum kinetic equation (0f|) 
for Eq — 3 x I0 8 V/cm and an initial temperature 
T = 20, 000K. Fig. g.a shows the result for a density 
of n = 10 22 cm -3 at different times (obtained by Fourier 
transforming j(i) over time intervals of 5 laser periods). 
One clearly sees the formation of a strong third harmon- 
ics which is comparable to and may even exceed the fun- 
damental component, i.e. 773 ~ 1. At later times, this 
peak vanishes again. We mention that higher odd har- 
monics are found for higher field strength also, besides, 
there is a clear signature of the second harmonics, cf. 
Fig. Ej. The temporal behavior of the harmonics emis- 
sion can be understood from the time evolution of the 
nonequilibrium cycle-averaged collision frequency 



v ei (E ,n;n, t) = 47T— 



n 2 j(t)E(t) 



(6) 



-pi E2(t) 

which is shown in Fig. ||. Due to the laser heating of the 
plasma, cf. Fig. |l|, and relaxation towards a Maxwellian, 
the collision frequency v e i decreases in time, limiting the 
efficiency T)pf(t) and, thus, the duration of the harmonics 
pulse. 
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FIG. 3. Nonequilibrium collision frequency, Eq. (B), at the 
initial moment, after 9fs and after 25 fs (white, grey and black 
symbols, respectively), n — 10 22 cm -3 , v t h = \/ {vi}- Due to 
collisional energy absorption, cf. Fig. 1, during the relaxation, 
symbols move to the left. Symbols with same shape refer to 
the same field strength, lines are guides to the eye. 

Finally, to understand the reason for the anomalously 
high efficiency of harmonics generation, we consider in 
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Fig. y the temporal evolution of the electron distribu- 
tion. The three figure parts show snapshots of / e (fc x , k z ) 
after 0,6 and 12 complete field cycles. Clearly, the dis- 
tribution is becoming anisotropic already during the first 
few laser periods. It strongly deviates from a Maxwellian, 
developing side peaks along the field direction, the dis- 
tance of which from the origin is proportional to the 
field strength. This peculiar nonequilibrium behavior is 
caused by electron-ion collisions in the presence of the 
field. It is observed for all considered densities and field 
strengths (see above), but is strongest around fl = u> p i 
and vq = Vth- Under these conditions, an electron spends 
the maximum possible time in the field of an ion [ex- 
tending over approximately ru = k^ 1 ] per laser cycle, 
^Tcoii ~ Qrn/vth ~ £l/u>pi ~ 1. These are the optimal 
conditions for collisional e-i momentum exchange in the 
field, as for electron acceleration (inverse BS), cf. Fig. [j], 
and harmonics emission (BS), Fig. ||.b. 

Our calculations show that, as a result of e-e collisions, 
the side peaks of f e start to smear out after 30 . . . 50fs, 
and after t le \ ~ 50 . . . lOOfs (depending on the density) a 
monotonically decaying distribution with a temperature 
of the order of several 10 e K is reached. By this time 
the plasma has become essentially collisionless: v e i has 
decayed to its equilibrium value flq] and v ee « 0, mak- 
ing it necessary to include collective and hydrodynamice 
(spatial heat flow, expansion etc.) effects, as mentioned 
in the Introduction, while for the treatment of collisions 
simpler models may be used. 



fs 10 fs 20 fs 




FIG. 4. Evolution of the electron distribution function in 
momentum space (the field is in k z direction). Figure shows 
snapshots at t = (left), after 6 laser periods (middle) and 
after 12 laser periods (right), when the main peak of f e is at 
the origin (the distribution as a whole oscillates with the field 
in fc z -direction and is isotropic in the fc^-fc^-plane). Note the 
different vertical scales. 

In summary, we have presented a selfconsistent nu- 
merical investigation of the early relaxation stage, t < 
50fs, of a pre-ionized dense plasma in a laser field with 
Eq = 10 7 . . . 10 9 V/cm.. We have found the formation of a 
strongly nonequilibrium electron distribution giving rise 
to coherent bremsstrahlung emission of ultrashort third 
harmonics pulses. An anomalously high conversion ef- 
ficiency (more than 6 orders of magnitude compared to 
plasmas in equilibrium) is predicted, making this funda- 
mental process observable in experiments. The most fa- 



vorable conditions are expected for a pre-ionized weakly 
undercritical gas plasma with high charge state (y e i ~ Z) 
excited by a linearly polarized laser pulse with a duration 
t p shorter than 50fs where the third harmonics should be 
observable in transmission. To further improve the quan- 
titative predictions, the plasma generation process has to 
be included in the simulation. For example, ionization 
from excited atomic levels should enhance the harmonics 
yield 0, accelerate the plasma heating and thus further 
reduce the duration of the higher harmonics pulse. Fi- 
nally, we mention that our results are not limited to laser 
plasmas, but are equally important for electron-hole plas- 
mas subject to intense THz radiation. 
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